Importance of convection in the compaction mechanisms of anisotropic granular media 
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We report the experimental observation of novel vortex patterns in a vertically tapped granular 
media. Depending on the tapping acceleration two behaviors are observed. For high acceleration 
a convection vortex appears in the whole media whereas for low acceleration two unstable vortices 
appear in the upper part of the media and slowly compact the lower part. We explain the formation 
of the vortices and relate them to granular convection. Our results demonstrate the importance of 
compression waves propagation on granular compaction. 

PACS numbers: 45.70.Cc, 45.70.Mg, 83.80.Fg 



5h 
43 



c 

o 
o 



> ■ 

00 ' 

: 

(N . 

o ■ 
in : 

o . 

-i— > ■ 
a ■ 



c 

o 
o 



X 



In hydrodynamics, flow instabilities lead to the forma- 
tion of vortex patterns which affect the behavior of the 
fluid. It is well known that granular media can exhibit 
both fluid-like and a solid-like behavior and thus do not 
behave like classical fluids. In order to understand the 
specificities of granular media, some studies of instabili- 
ties have been carried out for vertical vibratedgranular 
layers 0,0,0; granular flows and gazes @|. More- 
over, a granular packing under vertical shaking leads to a 
global compaction If tllflllll and this compaction can be 
linked to convection jlOl Ilia , Il3| . Since convection in flu- 
ids can lead to instabilities, can granular systems under 
compaction exhibit vortices ? In this letter we present 
a vortex creation process that occurs in a tapped pile of 
anisotropic grains and explain the mechanisms of such 
a creation. Since it has been recently proved that grain 
anisotropy is important in pattern formation Q we focus 
on anisotropic granular materials. 

The experimental setup consists of a glass cylinder of 
diameter D w 10 cm filled with 600 g of grains (corre- 
sponding to a height of roughly 10 cm). This container is 
placed on a plate connected to an electromagnetic exciter 
(LDS V406) which induces a vertical displacement of the 
plate. The container is, in this way, regularly shaken 
(At = Is) by vertical taps. Each tap is created by one 
entire period of sine wave at a fixed frequency / = 30 Hz. 
The resulting motion of the whole system, monitored by 
an accelerometer at the bottom of the container, is how- 
ever more complicated than a simple sine wave. At first, 
the system undergoes a positive acceleration followed by 
a negative peak with a minimum, —'y max . After the ap- 
plied voltage stops, the system relaxes to its normal re- 
pose position. This negative peak acceleration is used to 
characterize the tap intensity by the dimensionless accel- 
eration r = jmax/g (with g = 9.81 m.s~ 2 ). It should be 
pointed out that we only used taps for which the grains 
took off from the bottom of the glass cylinder i.e. above 
the lift-off threshold. The packing fraction is measured 
using a 7-ray absorption setup 01 • The results presented 
here have been obtained using long rice (basmati rice). 
In order to quantify the anisotropy, each grain is assim- 
ilated to an ellipsoid. The mean ratio between the two 
axes is 2.5 with 10% of dispersion. The length of the 



longer axe will be noted d. The method used to build 
our initial packing, already used for spheres |10| . is re- 
producible (we obtain the same initial packing fraction 
$ ~ 0.56 and the same packing fraction evolution for a 
given T). Sequences of 10 4 to 10 6 taps are carried out. 
The only control parameter is the tapping intensity T. 
Here we call " time" the number of taps and the " dynam- 
ics" is the succession of static equilibrium induced by the 
taps. The mean packing fraction increases monotically 
with time following stretched exponential. All the details 
on the compaction dynamics can be found elsewhere |12| . 

The first observation is that convection takes place in 
the medium : after r app , a given time depending on T, 
one or two counter-rotating vertical convection rolls ap- 
pear. They are easily identified by tracking grains posi- 
tion after each tap. To characterize them, we measure 
the time it takes a grain to revolve around the center 
of the roll, T conv . Fig. reports the variations of r conv 
and r app with T. As expected, the greater the tapping 
acceleration is, the smaller these two times are. Concern- 
ing the number of rolls, two cases should be considered. 
For high tapping intensity (typically T > 3) convection 
takes place within the whole medium. After about ten 
taps, two convection rolls appear but this situation be- 
comes unstable. One of the rolls progressively disappears 
and after this transient the whole medium is occupied by 
only one roll and the free surface of the medium is tilted 
from the horizontal (for example about 20° for T = 6). 
A sketch and a snapshot of this convection roll is re- 
ported FigUjD. The arrows represent in arbitrary units 
the displacement field of the grains and clearly demon- 
strate the existence of convection. A radial displacement 
profile is reported Fig^;. An immobile core can be ob- 
served. In an ideal system this non-mobile part should be 
reduced to a point, but due to the formation of spatially- 
correlated clusters the size obtained is around 4 d. After 
this part a linear profile is found with a slope compatible 
with the value of T conv . This part of the packing rotates 
around the immobile core like a solid. For lower values 
of the tapping intensity (typically T < 3) two convection 
rolls roughly equally-sized and localized near the free sur- 
face appear. The other part of this packing (below the 
rolls) is not submitted to convection. As for the high 
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FIG. 1: a) T app and T conv (see text for definition) versus 
tapping intensity V. b) Sketch of the convection obtained 
for basmati rice and for F = 6 (left) and snapshot of the 
medium during relaxation (right) where convection roll can 
be clearly seen on the walls, c) Radial displacement profile : 
the granular medium rotate like a solid around an immobile 



tapping intensity case, this situation is unstable. Con- 
trary to the previous case the two rolls do not merge : 
their size continuously decreases until they totally dis- 
appear. Note that the convection turns out not to be 
a phenomenon in competition with compaction. On the 
contrary the convection rolls create at their base a band 
of high packing fraction. On Fig. taken for T = 2.4 
after 1500 taps, this band is clearly visible at the base 
of the two convection rolls and one can observe that the 
grains in this band have a preferential horizontal orien- 
tation, due to the motion of the rice in convection. The 
band of high density can be distinguished on the packing 
fraction profile (FigEt)- The convection rolls create the 
high density zone and their size decrease corresponds to 
the growth of this zone. The orientation of the grains is 
totally different from the one observed in another com- 
paction experiment [T^j . Indeed, the vessel used in 0] 
is very narrow (1 grain length) compare to our (about 20 
grain lengths). Thus the nematic ordering observed in 
this tube is likely due to the strong confinement and no 
convection is possible ||. Note that the aspect ratio of 
the grains is probably also an important parameter. 
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FIG. 2: a) Sketch of the convection obtained for basmati rice 
and for V = 2.4 (left) and snapshot of the medium during 
relaxation (right). The compact zone can be observed on the 
side walls. b)Packing fraction profile obtained for basmati 
rice at F = 2.4 for various number of taps (see the arrow) : 
N = 0, 10, 100, 3116, 10000, 31162 and 100000. 



To understand the convection creation and evolution 
we have used a high-speed camera (1000 images per sec- 
ond) to analyze the motion and the deformation of the 
packing during a tap. Four phases can be visually dis- 
tinguished (see Fig. [3J : 

• the rice follows the motion of the plate (phase 1) 

• the rice takes off from the bottom (phase 2) 

• the rice lands on the bottom (phase 3) 

• a compression wave propagates throughout the 
medium (phase 4) 
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FIG. 3: Sketch of the four phases of a tap. The dashed line 
corresponds to the repose position of the vessel. The medium 
follow the motion of the vessel (phase 1), takes off from the 
bottom (phase 2), lands on the bottom (phase 3) and finally 
the compression wave propagates (phase 4). 

The rate of compaction and dilatation is measured dur- 
ing each phase for the first tap (the motion of the grains 
is greater here than during the other solicitations and 
thus the detection of motion is easier. The data were 
averaged on two different experiments and no significant 
differences were observed. To track the grains, we use an 
image processing software that computes the gray level 
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profile along a vertical line at the side wall. For clarity 
each profile is reported in the reference frame of the ves- 
sel. To improve contrast and thus to have sharp peaks on 
the gray level profile, 8% of the grains are paint in black. 
We can then follow the displacement of the peaks be- 
tween each picture (1000 pictures per seconds) and thus 
detect dilatation, compaction or wave propagation during 
each phase of the motion. The precision on the peak dis- 
placement is 0.3 mm for a packing height of about 10 cm. 
Note that, in order to avoid overlapping, each profile pre- 
sented here is shifted vertically. On FigQJi, the two gray 
level profiles correspond respectively to the state before 
and after the first tap of T = 6. These two profiles are 
similar and the peaks present can be easily tracked. Let 
us now examine with this method the phases of the mo- 
tion. During the first phase, the rice follows the motion 
of the container. As the medium is in a stable configura- 
tion, this motion does not create a relative displacement 
of the grains. As the cylinder slows down, the medium 
takes off from the bottom and a dilatation of the piling 
is observed (Fig^J)). This dilatation is not uniform : 1 
mm for the upper part of the packing and 0.5 mm for the 
lower part. A study of the orientation of the grains dur- 
ing this phase shows that they have a tendency to align 
vertically. Moreover, as already observed for glass beads 
in [F| , all the packing takes off and its bottom is still par- 
allel to the plate. During the third phase (Fig Et), the 
rice lands, which induces a compaction of the medium. 
Surprisingly this compaction is smaller than the dilata- 
tion monitored during the second phase and at this stage 
the medium is less compact than initially. Note that 
the vertical orientation tendency of the rice always exists 
even if it is less visible than during the second phase. In 
fact, the medium mainly compacts during the last phase 
(propagation of the compression wave) . The evolution of 
the gray level during that phase is reported on Fig 0Ji 
and compaction of the medium is indeed observed. The 
compression wave can be easily tracked on the fast cam- 
era movie. Its velocity can then be measured. For T = 6 
we have V ~ 16.5 m.s -1 . Moreover, it is only during this 
phase of propagation that the contribution to convection 
appears. Indeed, during phases 1, 2 and 3 the bottom 
of the packing remains parallel to the plate, hence no 
convection-like grain motions are observed. This result 
clearly demonstrates that compaction and convection are 
strongly linked. The compression wave is a consequence 
of the dissipation of the interaction energy stored by the 
grains during the previous motion phases. This allows 
the grains to move slightly around their initial position 
and to possibly find a more stable position corresponding 
to an higher value of the packing fraction. The velocity 
of the compression wave is of the same order of mag- 
nitude than those obtain in the same setup with glass 
beads [13J but smaller than the values obtained in a one 
dimensional granular medium 

It should be pointed out that in our set-up compaction, 
is mainly associated to compression wave propagation 



and convection. Nevertheless it is not the only one, 
since compaction can be observed even without convec- 
tion This also explain why the characteristic times 
found in are smaller than those obtained in 0. 
For r = 6, the effect of the compression wave is easy 
to analyze for two reasons. First, the dilatation of the 
medium just before the wave is still important and its 
effect is easier to appreciate. Second, the intensity of the 
compression wave is high because the medium lands on 
the plate when the plate comes up again. 
We have also performed such an analysis for V — 2.4. On 
Fig [SJi, the two gray level profiles reported correspond 
respectively to the state before and after the first tap. 
The motion of the grains is smaller than for T = 6. Many 
characteristics of the medium motion are similar to those 
obtained for T = 6. During the first phase (the grains fol- 
low the plate motion) the medium is not affected. Dur- 
ing the second phase, we monitor a small dilatation in 
the medium. Nevertheless the mean orientation does not 
change as much as for T = 6. During the third phase, the 
packing fraction increases. What is surprising and differ- 
ent from the r = 6 case is that during the last phase, the 
compression wave propagates inhomogeneously through 
the medium. Indeed no motion is visible at the bottom 
whereas a small displacement of 1 mm is monitored on 
the top of the packing. This is in agreement with previ- 
ous analysis : the wave does not affect the bottom part of 
the packing where no convection is observed, contrary to 
the upper part where both convection and compression 
wave propagation are present. Note that, as expected, 
the effect of the compression wave is less important for 
T = 2.4 than for T = 6. 

We explain this behavior difference by the motion of the 
packing relatively to the bottom plate. For T = 2.4, the 
medium lands as the plate is moving downwards. There- 
fore, the apparent fall velocity of the medium is smaller 
than in motionless case and the intensity of the compres- 
sion wave is smaller than for r = 6. The other observa- 
tion is that there is also a strong correlation between the 
apparition of convection and of the compression wave. In- 
deed, in all our experiments the compression wave creates 
a motion only in the part of the medium where convec- 
tion will take place during the other solicitations. This 
explanation helps understanding the transition between 
the two behaviors : for T = 3, the fall of the medium co- 
incides with the motionless, lowest position of the plate. 

To summarize, the new instability presented in this 
paper leads to a vortices formation during compaction 
of anisotropic granular media at relatively low tapping 
intensity. Below the vortices, an ordered compact zone 
can be created where the grains are oriented horizon- 
tally. By analyzing the packing behavior during one tap 
we showed that the main compaction mechanism is due 
to a compression wave propagation. Furthermore this 
compression wave propagation creates convection prov- 
ing the strong correlation between granular compaction 
and convection. 
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FIG. 4: Gray level profiles along a vertical line on the picture of the system for F — 6 : a) before (down) and after (up) the first 
tap, b) just before the medium takes off (down) and at its higher position, c) at its higher position (down) and just after the 
medium lands (up), d) just before (down) and after the compression wave. For clarity, each curve has an offset in the y-axis 
and is reported in the reference frame of the vessel. 
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FIG. 5: Gray level profiles along a vertical line on the picture of the system for F = 2.4 : a) before (down) and after the first 
tap (up), b) just before the medium take off (down) and at its higher position (up), c) at its higher position (down) and just 
after the medium lands (up), d) just before (down) and after (up) the compression wave. For clarity, each curve has an offset 
in the y-axis and is reported in the reference frame of the vessel. 



